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A problem of the equivalence of statistical ensembles is critically analyzed. It is shown, that 
although different probability distributions of statistical physics have the same behavior in the 
thermodynamic limit, there are physical observables - semi-intensive variables - which keep memory 
of the underlying ensembles. This property is an universal one and can be observed even in the 
simplest case of the grand canonical and canonical ensembles of the classical statistical physics. 



o 

^ 1 I. INTRODUCTION 

There is a standard textbook knowledge about the equivalence of different statistical ensembles in the 

■ thermodynamic limit. It is usually understood that in the thermodynamic limit averages of intensive variables 
calculated with different ensembles have the same values. It's also seemed an obvious corollary that one gets 
the same limiting expressions of all relevant physical quantities. So it was quite surprising result [ij that 
scaled particle variance variables defined as 

(N) 1 ' 

with N being the number of particles of a given charge, are different in the grand canonical, canonical 
and microcanonical ensembles - even in the the thermodynamic limit. Corresponding ensembles considered 
there were canonical or grand canonical with respect to the conserved charge, while microcanonical means 
here an exact energy and charge value of the system. Those results were also repeated for a more realistic 
hadron- resonance gas Q . 

■ The statistical equivalence problem is not of academic interest only. Thermal statistical models are 
| simple and effective tool to describe particle production in high energy heavy ion collision Q. A given 

t-H statistical ensemble is not a matter of an arbitrary choice but is related to internal properties of the system 
under consideration. Different statistical ensembles for finite size systems lead to easily detected different 
. predictions. However ,in the case of ultrarelativistic heavy ion collision departures from the thermodynamic 
<r— ( ' limit can be ignored [4| . 

Results of [J triggered discussions about statistical equivalence problem. It was however rigorously 
shown Q that appropriate probability distributions have the same functional limit in the thermodynamic 
limit. This explain the notion of the statistical ensemble equivalence problem. It has also appeared Q 
that there is a general class of physical quantities which are finite in the thermodynamic limit but are still 
different for different statistical ensembles. Those quantities were named semi-intensive variables. The scaled 
variance |T]) which was considered in [J, 0] belongs also to that class. 

Let us consider here the simplest example of the standard statistical physics. It will be shown that even 
in this case the notion of semi-intensive quantities is relevant for the physical situation. 



II. CANONICAL AND GRAND CANONICAL ENSEMBLES OF STATISTICAL PHYSICS 

The canonical ensemble is characterized by the fixed number of particles N, so the corresponding proba- 
bility distribution is just Dirac delta function 

P%(N,V) = S(N - N) . (2) 

The grand canonical ensemble is characterized by the fixed averaged number of particles (N) . The number 
of particles is now a random variable. The probability distribution to have N particles in the system is 
denoted as V?£(N, V). 
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The thermodynamic limit is defined as the infinite volume limit taken under the condition that corre- 
sponding densities are kept constant. This means 

N 

V — » oo, TV — > oo, — = n , (3) 

for the canonical ensemble and 

V - oo, (iV) - oo, M = („) , (4) 

for the grand canonical ensemble. 

To formulate correctly the thermodynamic limit of quantities involving densities, one should define the 
probability distributions for densities. They are 

V%(n,V) = W$ n {Vn,V)=6{n-n), (5) 

for the canonical ensemble and 

V^(n,V)=VF^ n) (Vn,V), (6) 

for the grand canonical ensemble. 

The equivalence of the grand canonical and canonical distributions means 

7>&f°(n) = Hm Vffi(n, V) = 7#(n) . (7) 

In the thermodynamic limit the only relevant probabilities distributions are those related to densities. These 
distributions are expressed by probability moments. They are calculated for densities but not for particles. In 
the practice, however, we measure particles — not densities as we do not know related volumes. Fortunately, 
volumes can be omitted by taking corresponding ratios. 

Let us consider e.g. the density variance An 2 . This can be written as 

An 2 = (n 2 ) - (n) 2 = ■ ' v ' 



V 2 



By taking the relative variance 

An 2 (N 2 ) - (iV) 2 



(n) 2 (N) 2 

volume-dependence vanishes. 

A special care should be taken for calculations of ratios of particles moments. Although moments are 
extensive variables their ratios can be finite in the thermodynamic limit. These ratios are examples of semi- 
intensive variables. They are finite in the thermodynamic limit but those limits depend on volume terms 
in density probability distributions. One can say that semi-intensive variables "keep memory" where the 
thermodynamic limit is realized from. 

Let consider as an example the scaled particle variance 

(N 2 ) - (AO 2 = y < 



(N) 
The term 



tends to zero in the thermodynamic limit as ©(V -1 ). So a behavior of the scaled particle variance depends 
on the 0(V~ l ) term in the scaled density variance. A more detailed analysis of semi-intensive variables is 
given in Q. 
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III. POISSON DISTRIBUTION IN THE THERMODYNAMIC LIMIT 



To clarify this approach let us consider a well known classical problem of Poisson distribution but taken 
in the thermodynamic limit. 

Let us consider the grand canonical ensemble of noninteracting gas. A corresponding statistical operator 

is 



,-/3H+-iN 



D 



This leads to the partition function 

where z is one-particle partition function 

z(T,V) 



Z{V,T,i)= e ze \ 



V 



(2tt) 2 



d 3 p e~P E{p) = Vz (T) , 



(8) 
(9) 
(10) 



A 7 parameter (= f3/i) is such to provide the given value of the average particle number (N) = V(n). This 
means that 



=7 _ 



(11) 



The resulting probability distribution to obtain N particles under condition that the average number of 
particles is (N) is equal to Poisson distribution 



U (AT\ W_ -<N) 

We introduce corresponding probability distribution V for the particle number density n = N/V 

(V{n)) v * 



VOC {n . v) = VPv{n){Vn) = v S^l^ 



-V(n) 



The large volume behavior can be obtained here by means of the saddle point method. This gives 

Vn 



v ' V27rn V n 



,V(n-{n)) 



1 



1 



12Vn 



oiy- 



(12) 



(13) 



(14) 



This expression in singular in the V — » oo limit. To estimate a large volume behavior of the probability 
distribution (I13|) one should take into account a generalized function limit. So we are going to calculate an 
expression 



(G)v= I dnG(n)Vgf(n;V), 

for an arbitrary function G(n) and with V^{n\ V) replaced by the asymptotic form from Eq (fT4] 
An asymptotic expansion of this integral is given by the classical Watson-Laplace theorem 

Theorem Let I = [a, b] be the finite interval such that 

1. max5(s) is reached in the single point x — xq, a < xq < b. 

2. f(x),S(x)eC(I). 

3. f(x),S(x) e C°° in the vicinity of x , and S" (x ) ^ 0. 
Then, for A — > oo, A € 5 e , there is an asymptotic expansion 



F[X] ~ e xs ^J2 c * X ~ k ~ 1/2 

k=0 

r(fc + 1/2) / d 



2 k 



S(x ) - S(x )Y k - 1/2 
(x - x ) 2 



(15a) 
(15b) 



5 £ is here a segment | argz| ^ 5 — e < -| in the complex z-plane. 
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In the next to leading order asymptotic expansion one gets then 

(G) v = G((n)) + { ^G"({n)) + 0{V- 2 ) , (16) 

for any function G. 

This gives us the exact expression for the density distribution (|13| in the large volume limit 

Vffifa V) = S(n - (n» + M S "(n - (n» + 0{V~ 2 ) . (17) 

This confirms the equivalence of the grand canonical and canonical distributions as was defined in Eq ([?]). 
We are now able to obtain arbitrary density moments up to 0(V~ 2 ) terms. 

,fc\ _ f j tn /„.T/\_/„\fc i - 1) / \fe-i i /nnr-2\ 



(n k ) v = I dnn k V {n) (n-V) = (n) k + v — ' (n)^ 1 + 0(V~ 2 ) . (18) 



We have for the second moment (intensive variable!) 



This means 



An 2 = ^^0. (19) 

as expected in the thermodynamic limit. 

The particle number and its density are fixed in the canonical ensemble so corresponding variances are 
always equal to zero. The result (I19p can be seen as an example of the equivalence of the canonical and 
grand canonical distribution in the thermodynamic limit. This equivalence is clearly visible from the Eq 
(I17p where the delta function in the first term can be considered as the particle number density distribution 
in the canonical ensemble. 

A more involved situation appears for particle number moments (extensive variable!). Eq (|18[) translated 
to the particle number gives 

(N k ) = (N) k + k<yk - ^ (N) k ~ 1 + 0(V k - 2 ) , (20) 



One gets for the scaled variance (semi-intensive variable!) 

AN 2 



(AO 



(21) 



what should be compared with zero obtained for the canonical distribution. 

The mechanism for such a seemingly unexpected behavior is quite obvious. The grand canonical and 
the canonical density probability distributions tend to the same thermodynamic limit. There are different 
however for any finite volume. Semi-intensive variables depend on coefficients at those finite volume terms 
so they are different also in the thermodynamic limit. 

Similar calculations can also be done for the energy taken in the canonical and in the grand canonical 
distribution (7|. A conclusion is also similar. Energy probability distributions coincide in the thermodynamic 
limit for both ensembles, although finite volume corrections are different. 

A. Semi-intensive variables 

Now we are in position to create some semi-intensive variables. They are finite in the thermodynamic limit 
and have different values dependently on how the charge conservation is implemented in the description of 
the system. There is actually a broad class of variables. We take an an example 

(N k ) (N) k 

Sk = (jv)*-i ■ ( } 
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Taking into account Eq ([20]) one gets in the thermodynamic limit 

The scaled variance is just a special case of corresponding to fc = 2. 

Another examples are quantities closely related to cumulant or factorial cumulant moments or suscepti- 
bility ratios. Let define p— th order susceptibility 

_ d p \nZ 

where Z is a general partition function - not only the Poissonian one. 
One can easily check that the ratios 

fcp;r — ) (^) 

are semi-intensive quantities. 

One can also construct more involved semi-inclusive variables having a finite thermodynamic limit behavior 
which are determined by higher order asymptotic terms of the corresponding probability distributions. 



IV. CONCLUSIONS 



We have discussed the differences in the asymptotic properties of the probability functions for a very 
simple classical system with both exact and average particle number conservation. We have shown that in 
the thermodynamic limit the corresponding probability distributions in the grand canonical and canonical 
ensembles coincide and are described as generalized functions. This property is a direct consequence of the 
grand canonical and canonical ensemble equivalence in the thermodynamic limit. However, the first and 
higher finite volume corrections to the asymptotic value differ in these ensembles. 

Finally, we have derived the asymptotic behavior of the particle moments using the results obtained for 
the probability functions. We have also applied these results to find the thermodynamic of semi-intensive 
quantities. It was shown that in systems with exact and average particle number conservation such quantities 
should naturally converge to different values in the thermodynamic limit. This is because the behavior of the 
semi-intensive quantities in the near vicinity to the thermodynamic limit are determined by the subleading, 
finite volume corrections to the probability distributions which are specific to a given statistical ensemble. 

Semi-intensive variables are well suited to become benchmarks of different statistical distributions - 
even in the thermodynamic limit. That property is an universal one and, as we've seen, it appears also 
on the simplest level. Is should be noted here that first moments of basic physical quantities tend to the 
same thermodynamic limit for different statistical ensembles. Equations of state are constructed as relations 
between first moments so they do not provide any information about underlying statistical ensembles. It can 
be achieved only by semi-intensive variables as e.g. scaled variance or other fluctuations related quantities. 
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